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Electroweak Baryogenesis and the Higgs and Stop masses * 

Mariano Quir6s a 

a Instituto de Estructura de la Materia, Serrano 123 E-28006, Spain 

In this talk we review the actual situation concerning electroweak phase transition and baryogenesis in the 
minimal supersymmetric extension of the Standard Model. A strong enough phase transition requires light Higgs 
and stop eigenstates. For a Higgs mass in the range 110-115 GeV, there is a stop window in the range 105-165 
GeV. If the Higgs is heavier than 115 GeV, stronger constrains are imposed on the space of supersymmetric 
parameters. A baryon-to-entropy ratio is generated by the chargino sector provided that the fi parameter has a 
CP-violating phase larger than ~ 0.04. 



1. Introductory remarks 

Electroweak baryogenesis [Q is an appealing 
mechanism to explain the observed, Big Bang 
Nucleosynthesis (BBN), value of the baryon-to- 
entropy ratio 0, 77BBN = tib/s ~ 4 x 10~ 10 , at 
the electroweak phase transition ||, that can be 
tested at present and future high-energy collid- 
ers. Although the Standard Model (SM) contains 
all the necessary ingredients |Q for a successful 
baryogenesis, it fails in providing enough baryon 
asymmetry. In particular it has been proven by 
perturbative [Q and non-perturbative || methods 
that, for Higgs masses allowed by present LEP 
bounds P|, the phase transition is too weakly 
first order, or does not exist at all, and any pre- 
viously generated baryon asymmetry would be 
washed out after the phase transition. On the 
other hand the amount of CP violation arising 
from the CKM phase is too small for generating 
the observed baryon asymmetry M. Therefore 
electroweak baryogenesis requires physics beyond 
the Standard Model at the weak scale. 

Among the possible extensions of the Stan- 
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dard Model at the weak scale, its minimal super- 
symmetric extension (MSSM) is the best moti- 
vated one. It provides a technical solution to the 
hierarchy problem and has deep roots in more 
fundamental theories unifying gravity with the 
rest of interactions. As for the strength of the 
phase transition |S, a region in the space of su- 
persymmetric parameters has been found fl-^Sj 
where the phase transition is strong enough to 
let sphaleron interactions go out of equilibrium 
after the phase transition and not erase the gen- 
erated baryon asymmetry. This region (the so- 
called light stop scenario) provides values of the 
lightest Higgs and stop eigenstate masses which 
are being covered at LEP and Tevatron colliders. 

The MSSM has new violating phases that can 
drive enough amount of baryon asymmetry. Sev- 
eral computations have been performed p6[-[p8| 
in recent years, showing that if the CP-violating 
phases associated with the chargino mass param- 
eters are not too small, these sources may lead to 
acceptable values of the baryon asymmetry. In 
this talk I will present some aspects of a recent 
computation |39) of the CP-violating sources in 
the chargino sector which improves the computa- 
tion of Ref. in two main aspects. On the one 
hand, instead of computing the temporal compo- 
nent of the current in the lowest order of Higgs 
background insertions, we compute all current 
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components by performing a resummation of the 
Higgs background insertion contributions to all 
order in perturbation theory. The resummation 
is essential since it leads to a proper regulariza- 
tion of the resonant contribution to the temporal 
component of the current found in Ref. |27|] and 
leads to contributions which are not suppressed 
for large values of the charged Higgs mass. On 
the other hand we consider, in the diffusion equa- 
tions, the contribution of Higgsino number violat- 
ing interaction rate from the Higgsino fi term 
in the lagrangian, T^, that was considered in our 
previous calculations in the limit T^/T — > oo. 



Table 1 

Mass spectrum along the U direction. 



field 


d.o.f. 


mass 2 


4 gluons 


12 


9 2 U 2 /2 


1 gluon 


3 


2g 2 U 2 /3 


1 B gauge boson 


3 


g' 2 U 2 /9 


5 squark-goldstones 


5 


m 2 v +g 2 C/ 2 /3 


1 squark 


1 


mf r + g 2 3 U 2 


4 Q L -Higgs ^ 


4 


~m 2 H /2 + h 2 sin 2 /3C/ 2 


2 Dirac fermions (tj, , H ) 


8 


M 2 + h 2 U 2 



2. The phase transition 

The possibility of achieving, in the MSSM, a 
strong-enough phase transition for not washing 
out any previously generated baryon asymmetry 
characterized by the condition 

v(T c )/T c >l , (1) 

has been recently strengthened by three facts: 

• The presence of light t r (with small mixing 
At) considerably enhances the strength of 
the phase transition ||||]. This is the so- 
called light stop scenario. 

• Two-loop corrections enhance the phase 
transition in the SM 0, and in the 
MSSM 

• The validity of perturbation theory, and in 
particular the results of two-loop calcula- 
tions, for the light stop scenario has been 
recently confirmed by non-perturbative re- 
sults In fact, non-perturbative re- 
sults lead to a stronger phase transition 
than perturbative ones, by ~ 10 — 15%. 

However the price the light stop scenario has 
to pay is that it may require moderately neg- 
ative values of the supersymmetric parameter 
my = —friij and then, apart from the electroweak 
minimum along the Higgs {<p) direction, another 
color breaking minimum along the U = tR di- 
rection might appear. Therefore both directions 
should be studied at finite temperature. 



To systematically analyze the different possi- 
bilities we have computed the two-loop effective 
potential along the <f> and U directions and com- 
pared their cosmological evolutions with T. The 
two-loop effective potential along the </> direction 
was carefully studied in Ref. Jig] . The one-loop 
correction is dominated by the exchange of the 
top/stop sector while the two-loop effective po- 
tential is given by two-loop diagrams with stops 
and gluons, as well as one-loop diagrams with the 
stop thermal countcrterm. The two-loop effec- 
tive potential along the U direction was studied 
in Refs. p7| . The mass spectrum is given in Table 
1, where we have considered small mixing A t /m,Q 
and large gluino masses. The one-loop diagrams 
correspond to the propagation of gluons, squarks 
and Higgses as well as Dirac fermions. The lead- 
ing two-loop contributions correspond to sunset 
and figure-eight diagrams with the fields of Table 
1 propagating, as well as one-loop diagrams with 
thermal counterterm insertions corresponding to 
gluons and squarks. 

For a given value of the supersymmetric param- 
eters we have computed T c , the critical temper- 
ature along the (j> direction, and , the critical 
temperature along the U direction. Therefore, 
four different situations can arise: 

• a) Tj 7 < T c and rhu < = 

(mf I v 2 g^/12) 1 ^ 4 |J. In this case the phase 
transition proceeds first along the (f> direc- 
tion and the field remains at the electroweak 
minimum forever. This region is called sta- 
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bility region. 

• b) < T c and m,jj > mfj. In this region 
the electroweak minimum is metastable 
(metastability region). It can be physi- 
cally acceptable provided that its lifetime 
is larger than the age of the universe at this 
temperature: < H. 

• c) T c c/ > T c and fhu < m u . In this case the 
[/-phase transition happens first and there- 
fore it would be physically acceptable pro- 
vided that the lifetime is shorter than the 
age of the universe: IV^ > H. This re- 
gion is called two step region and has been 
excluded [|||. 

• d) T c y > T c and fhu > mfj. This is the 
region of instability of the electroweak min- 
imum. It is absolutely excluded. 




m H (GeV) 



Figure 1. The absolute region of stability in the 
(niH,mj R ) plane, for wlq — 2 TeV. The allowed 
region is inside the solid lines. 

Fig. [j] shows the absolute region of stability for 
rriQ < 2 TeV, where we can see an absolute upper 
bound on the Higgs mass ~117 GeV. In view of 
the non-perturbative results on the phase tran- 
sition H, we have used in Fig. Ill the condition 



v(T c ) Z 0.9T C in our perturbative two-loop calcu- 
lation. 

3. The CP- violating chargino currents 

Our aim in this section is to compute the Green 
functions^for the charged gaugino-Higgsino sys- 
tem h c -W c , describing the propagation of these 
fcrmions in the presence of a bubble wall. The 
bubble wall is assumed to be located at the space- 
time point z, where there is a non-trivial back- 
ground of the MSSM Higgs fields, H t (z), which 
carries dimensionful CP-violating couplings to 
charginos. We shall use these Green functions to 
compute the left-handed and right-handed cur- 
rents corresponding to 

«*>=(!;)• **>-(£)•<*> 

at the point z. 

Since the mass matrix depends on the space- 
time coordinates, and we must identify the free 
and perturbative parts out of it, in order to make 
such a selection we will expand the mass matrix 
around the point = (r, t) (the point where we 
are calculating the currents in the plasma frame) 
up to first order in derivatives as, 

M(x) = M(z) + (x- zYM^z) , (3) 

where we use the notation M^z) = dM(z)/dz ,J '. 
The chargino mass matrix in (ph is given by 

M(z)=( M > (4) 
\ ui{z) n c ) 

where we have defined Ui(z) = gHi(z). The mass 
eigenvalues are given by 

. (A + A + u 2 1 (z))M 2 + u 1 (z)u 2 {z)n* 

mi z = , 

V(A + A)(A + A) 

(A + A-«f(z)) fi c - ui(z)u 2 (z)M 2 

m 2 (z) — ■ . 

V(A + A)(A + A) 

where field redefinitions have been made in order 
to make the Higgs vacuum expectation values, as 
well as the weak gaugino mass M 2 , real, and 

A = (M 2 -\fi c \ 2 -u 2 +u 2 )/2 

A = (Ml -| Mc | 2 ~u 2 2 +u\)/2 

1 /2 

A = (A 2 + |M 2 ui+^> 2 | 2 ) " (5) 
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The vector and axial Higgsino currents can now 
be defined as: 

f Hh (z) - lim {Tr [P 2 ^S RR (x,y;z)] 

± Tr[P 2 a»S% L (x,y;z)]} (6) 

where P 2 = ((Jq — 03)/2 is a projection opera- 
tor and Sih{%, y, z) are the Green functions in the 
weak eigenstate basis after making the expansion 
(H) and a resummation to all orders in M(z). 

The detailed calculation of the currents can be 
found in Ref. |39f| . In particular only the time 
component is relevant for the diffusion equations 
of next section where we will use as sources of 
Higgs densities -y H ,h ^ r~ j% jh (z), and T~ ~ 
aw T is the inverse typical thermalization time. 
In this way we obtain for the vector and axial 
sources the expressions: 



1 f°° 
H(z) = —Re J dp (1 + 2/) 



1 



1 



Z\ Z 2 \Zl + Z 2 



(10) 



K(z) 



-^Re / dp (1 + 2/) 

^ JO 



1 



1 



z\ Zi \zx + z 2 



(11) 



/ = — np(|f> |), where np is the Fermi-Dirac dis- 
tribution function and z, is defined as 

Zi(p°) = y p o(p° + 2ir g ) -|mi(z)|2 (12) 



1H ~-2^ 5 2 r~Im(M 2Mc ) (7) 

{H\z)P'{z) [F{z)+Q{z)\ 
+g 2 H 2 (z) cos 2(3(z) [H(z)H'(z) sin 20(z) 
+ ff 2 (z)cos2/3(z)/3'(z)] W(z)} 

j h ~ 2w w5 2 r~Im(M 2Alc ) 
[i?(z)iJ'(z) sin2/3(z) + 7J 2 (z) cos 2/3(z)/3'(z)] 
{/C(z) + 2(A + A)W(z)} 

where the functions J 7 , Q, Ti, /C are defined as: 



^(z) = ^Ref oo dp°(l + 2f) 
1 ^ 3 



Zl + z 2 



(8) 



= ^Re jf dp //' 



1 



\zi+z 2 / |toi(z)| 2 - |m 2 (z)| 2 

Zi 

|toi(z)| 2 - |m 2 (z)| 2 -4ir~ p° 

h ^ 

|mi(z)| 2 — |m 2 (z)| 2 + 4ir~ p° 



(9) 




Figure 2. Plot of the sources jH,h, for the val- 
ues of supersymmetric parameters specified in the 
text, as functions of z/L u . 



with positive real and imaginary parts satisfying 
Re(zj) = T~ p°/Im(zi). Notice that the sources 
are proportional to the wall velocity v w , and so die 
when the latter goes to zero, which is a physical 
requirement. 

In Fig. H we plot the sources jH.h for a chosen 
set of supersymmetric and bubble parameters. In 
particular the wall velocity is chosen as v w — 0.05 
and the bubble wall width as L u = 20 /T which is 
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suggested by detailed numerical analyses of bub- 
ble formation [^9pc| l. For the supersymmetric 
parameters we choose uiq = 1.5 TeV, A t = 0.5 
TeV, M 2 = M = 200 GeV, tan/? = 20 and three 
different values of m A = 100, 150, 200 GeV, cor- 
responding to the dashed, solid and dot-dashed 
curves of Fig. ^. From it we can see two main 
features. On the one hand, jh is very sensitive 
to the value of itia, and so to the correspond- 
ing value of A/3, as expected, and decreases when 
via increases. On the other hand, 7^ is domi- 
nant only for large values of tua, such that the 
A/3 suppression of jh is stronger. But it is never 
overwhelming the contribution of 7^, because it 
is tan /3 suppressed. 

We also computed the bubble solutions of the 
MSSM |l{| using the two-loop effective potential 
and we have confirmed the goodness of the thick 
wall approximation (L U T C 3> 10 VtS). 




300 

m A (GeV) 



Figure 3. A/3, in units of 10 -3 , as a function 
of itia for the values of the supersymmetric and 
bubble parameters indicated in the text. 



In Fig. §, we plot A/3(T) for the following val- 
ues of the supersymmetric: M2 = fx = 200 GeV, 
A t = 500 GeV, m Q = 1.5 TeV. Notice that the 
value of T c , as well as A/3, is different from point 
to point, for different values of itia- We see from 



Fig. H that A/3 goes from ~ 10 2 for itia = 100 
GeV, to - 10~ 3 for m A = 500 GeV. This fact will 
have an important impact on the mechanism of 
baryon asymmetry. 

4. The diffusion equations 

To evaluate the baryon asymmetry generated 
in the broken phase we need to first compute 
the density of left-handed quarks and leptons, 
til, in front of the bubble wall (in the symmet- 
ric phase). These chiral densities are the ones 
that induce weak sphalerons to produce a net 
baryon number. Since, in the present scenario, 
there is essentially no lepton asymmetry, the den- 
sity to be computed in the symmetric phase is 
n L = n-Q + Y^i=i n Qi where the density of a chi- 
ral supermultiplet Q = (q, q) is understood as the 
sum of densities of particle components, assum- 
ing the supcrgauge interactions to be in thermal 
equilibrium, hq = n q + fiq. If the system is near 
thermal equilibrium, particle densities, n,i, are re- 
lated to the local chemical potential, fj,i by the 
relation rii — ki^iT 2 /6, where ki are statistical 
factors equal to 2 (1) for bosons (fermions). In 
fact, assuming that all quarks have nearly the 
same diffusion constant it turns out that | p6[ , 
n L = 5n Q + 4 T^T- 
One can write now a set of diffusion equa- 
tions involving tiq, tit, tlh 1 (the density of H\ = 

(h\, hi)) and nu 2 (the density of H2 = (^2,^2)), 
and the particle number changing rates and CP- 
violating source terms discussed above. In the 
bubble wall frame, and ignoring the curvature of 
the bubble wall, all quantities become functions 
of z = r + v u t. In the limit of fast Yukawa cou- 
pling Ty and strong sphaleron T ss rates, we can 
write the diffusion equations as: 



\n'q + 2 n' T — n' H ] = D q [h'q + 2 n?] 

n H 



r, 



v u [n' Q + 2n' T - n' h ] = D q \n!' Q + 2 n£] 
-D h < + r. 



TLT 
kx 



K H 



■ Ih 



(13) 



6 



where uq and tit are replaced by the (approxi- 
mate) explicit solutions 



tit = 



k Q (9k T - k B ) 
k H {k B + 9k Q + 9kr) 
k T (9k Q + 2k B ) 



k H (k B + 9k Q + 9k T ) 



(n H + 71,0(14) 



D q ~ 6/T and Dh ~ 110/T are the correspond- 
ing diffusion constants in the quark and Higgs 
sectors, n B = Uh 2 + n Bl , nh = tih 2 — n Hi , 
k B = fciTj + kn 2 , and r M corresponds to the 
fi c HiH2 term in the Lagrangian. There are also 
the Higgs number violating and axial top num- 
ber violation processes, induced by the Higgs self 
interactions and by top quark mass effects, with 
rates and r m , respectively, that are only active 
in the broken phase. 



, = 150 GeV 

r -o.i t 




symmetric phase (z < 0), where weak sphalerons 
are active, which is essential for the left-handed 
quark asymmetry to bias the weak sphalerons to 
violate baryon number. In particular, we can see 
from Fig. [| that the density n B is larger than 
n/i, for r M = 0.1 T. In previous analyses 26 27| 
the limit — > oo was implicitly assumed and nh 
was completely neglected. Fig. |] shows that this 
is not such a bad approximation. Furthermore, 
the relative importance of n B and is shown 
for different values of T M = 0.01 T, 0.1 T, T in 
Fig. | We can see that while for T p = 0.01 T 
rih is sizeable, for = T it is negligible in good 
agreement with our previous results p7[. 



m. = 150 GeV 




z/L 



z/L„ 



Figure 5. Plots of densities-to-entropy ratios 
n B ,h/s as functions of z/L u for different values 

ofr„, 



Figure 4. Plot of the different density-to-entropy 
ratios for the values of supersymmetric parame- 
ters specified in the text. 



The system of equations (|l3|) has been solved 
numerically in Ref . 39 1 , where also a good enough 
analytical approximation is provided. In Fig. |J we 
plot the numerical solution corresponding to the 
same set of supersymmetric and bubble parame- 
ters as in Fig. @. All densities diffuse along the 



5. The baryon asymmetry 

In this section we present the numerical results 
of the baryon asymmetry computed in the pre- 
vious section and, in particular, of the baryon- 
to-entropy ratio r\ = n B /s, where the entropy 
density is given by s — ^-g e ff T 3 , with g e ff be- 
ing the effective number of relativistic degrees of 
freedom. 
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Since we assume the sphalerons are inactive in- 
side the bubbles, the baryon density is constant in 
the broken phase and satisfies, in the symmetric 
phase, an equation where til acts as a source |26| 
and there is an explicit sphaleron-induced relax- 
ation term |p9| 

v LJ n' B (z) = -6(—z) [n F T ws n L (z) + Un B (z)] (15) 
where np = 3 is the number of families and 1Z = 



M.= n 



jnpT ws is the relaxation coefficient. Eq. (15) 



can be solved analytically and gives, in the broken 
phase z > 0, a constant baryon asymmetry, 



n B 



(16) 



The profiles H(z), (3{z) have been accurately 
computed in the literature [jl9| . For the sake of 
simplicity, in this work we use a kink approxima- 
tion 



H(z) 



1 — tanh 



p(z) = p - - A/3 [ I 



tanh 



i-i£ 



2z 
1 

L U1 



(17) 



This approximation has been checked to repro- 
duce the exact calculation of the Higgs profiles 
within a few percent accuracy, provided that we 
borrow from the exact calculation the values of 
the thickness L u: /2a and the variation of the an- 
gle P(z) along the bubble wall, A/3, as we will do. 
In particular we will take a = 3/2, L u = 20/T, 
and we have checked that the result varies only 
very slowly with those parameters, while we are 
taking the values of A/3 which are obtained from 
the two-loop effective potential used in our calcu- 
lation. 

In Fig. ^| we have fixed r/ = ??bbn and plot 
sin<^, where ip^ is defined as [i c = /xexp(i<p M ), 
as a function of /.i. We have fixed all bubbles 
and supersymmetric parameters as in Fig. g, fixed 
Mi = /-J. and ran over three typical values of 
the pseudoscalar Higgs mass tua- In all cases 
the phase transition is strong enough first order, 
v{T c )/T ~ 1, the running mass of the lightest 



- 100 GeV 

- 150 GeV 

- 200 GeV 




(GeV) 



Figure 6. Plot of sint^ as a function of M2 = 
H for the values of supersymmetric parameters 
specified in the text. 



stop is around 120 GeV and the Higgs mass is, 
within the accuracy of our calculations, between 
110 and 115 GeV. 

Since the computed rj behaves almost linearly 
in sint/? M , and we have fixed r\ = ijbbNi the more 
baryon asymmetry is generated the smaller the 
value of sint^. This can be seen from Fig. || 
where we have been working at the resonance 
peak Mi = /i and hence baryogenesis has been 
maximized. We see, for all values of uia, that the 
minimum of sin (p^ sits around /i ~ 200 GeV. The 
value of sin (p^ at the minimum decreases with 
tua- In particular, for uia — 100 GeV, which is 
about the lower limit from present LEP data (in 
fact, the present preliminary bounds from LEP |(| 
are tua > 89.9 GeV for large values of tan/3, as 
those we are using in our plots) we obtain that 
sin^ 0.04. The dependence of sin</? M with re- 
spect to supersymmetric and bubble parameters, 
as M2 ^ fJ-, vriA and has been thoroughly ana- 
lyzed in Ref. |3!| where the reader can find it as 
well as details of the calculation of supersymmet- 
ric chargino sources. 
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6. Conclusions 

The main conclusion after a detailed analysis 
of both the phase transition and the baryogenesis 
mechanism in the MSSM is that it is still alive af- 
ter the recent experimental results at high-energy 
colliders and, in particular, at the LEP collider. 

Concerning the phase transition, its strenght 
is controlled mainly by the Higgs mass and the 
lightest (right-handed) stop mass. Bubbles are 
formed with thick walls (the value of the thick- 
ness is ~ L u /3, with L u ~ 20/T) and prop- 
agate with extremely non-relativistic velocities 
(Vui ~ 0.1 — 0.01). The strenght of the phase 
transition has to be such that v(T) > T at the 
critical temperature. This imposes a strong con- 
straint on the supersymmetric parameters in or- 
der to avoid sphaleron erasure, in particular in 
view of the most recent (preliminary) bounds on 
the SM-like Higgs mass, tuh > 113.2 GeV and 
the observed excess of events with bb invariant 
mass ~ 114 GeV. Here two possibilities can be 
drawn: 

• The observed excess of events corresponds 
to a Higgs signal. 

In this case the combined Higgs mass and 
BAU requirements impose some restric- 
tions on the supersymmetric parameters. 

a) Heavy pseudoscalars and large tan/3: 
say vtia > 150 GeV and tan/3 > 5; 

b) Heavy left-handed stops and controlled 
stop mixing: tuq > 1 TeV and 0.25 5 
At/mQ ^ 0.4; and, c) Light right-handed 
stops: 105 GeV < m~ < 165 GeV. In this 
case the first prediction of the BAU scenario 
would have been realized and we would need 
confirmation for the rest, in particular for 
the light stop. 

• The observed excess of events does not cor- 
respond to a Higgs signal. 

In that case a reduction of the Hbb coupling 
would be needed. For the values of A t and 
\x consistent with electroweak baryogenesis, 
a reduction of the coupling of the CP-even 
Higgs boson to the bottom quark would de- 
mand not only small values of vtla — 100- 
150 GeV, but also large values of tan/3 > 10 



and of \fiAt\/rriQ > 0.1 (the larger tan/3, 
the easier suppressed values of the bottom 
quark coupling are obtained). A detailed 
discussion on this issue has already been 
done |§. 

Finally, concerning the generated baryon asym- 
metry, we have found that it requires the CP- 
violating phase to be, ip^ > 0.04. Values of 
ipp £ 0.04 can lead to acceptable phenomenol- 
ogy if either peculiar cancellations in the squark 
and slepton contributions to the neutron and elec- 
tron electric dipole moments (EDM) occur H), 
and/or if the first and second generation of 
squarks are heavy. This second possibility is 
quite appealing and, as has been recently demon- 
strated [ 42jl . leads to acceptable phenomenology. 
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